EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



Oh! 



> 

O 
(N 

O 

a^ 
o 
o 
o 



CERN-EP-2000-100 
July 20*"^ 2000 



O 

§ : A Measurement of .the Rate of Charm 

(n: Production in W Decays 



The OPAL Collaboration 



Abstract 



Using data recorded at centre-of-mass energies around 183 GeV and 189 GeV with the OPAL detector 

^ , at LEP, the fundamental coupling of the charm quark to the W boson has been studied. The ratio 

I \ Tx^ = r (W -^ cX)/r(W -^ hadrons) has been measured from jet properties, lifetime information, 

ij'' and leptons produced in charm decays. A value compatible with the Standard Model expectation of 

,^ . 0.5 is obtained: R^ = 0.481 ±0.042 (stat.) ±0.032 (syst.). By combining this result with measurements 

of the W boson total width and hadronic branching ratio, the magnitude of the CKM matrix element 

L^ ■ iV^sl is determined to be WcJ = 0.969 ± 0.058. 
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1 Introduction 



Since 1997 the LEP e'''e~ collider has been operated at energies above the threshold for W-pair 
production. This offers a unique opportunity to study the hadronic decays of W bosons in a clean 
environment and to investigate the coupling strength of W bosons to different quark flavours. The 
fraction of W bosons decaying hadronically to different quark flavours is proportional to the sum of 
the squared magnitudes of the corresponding elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix [|| . A measurement of the production rates of different flavours therefore gives access to the 
individual CKM matrix elements. 

In e'''e~ -^ W^W~ events, the only heavy quark commonly produced is the charm quark. The 
production of bottom quarks is in fact highly suppressed due to the small magnitude of l^ubl and | Vcbl 
and the large mass of the top quark, the weak partner of the bottom quark. This allows for a direct 
measurement of the production fraction of charm in W decays without a separation of charm and 
bottom quarks. The magnitude of the CKM matrix element V^s can then be derived from the charm 
production rate, using the knowledge of the other CKM matrix elements. So far direct measurements 
of |V^s| have limited precision compared to the other CKM matrix elements important in W decays, 
i.e. \Vad\, iKisI, and |V^d|- The most recent evaluation of the magnitude of V^g from exclusive charm 
semielectronic decays yields |V^s| = 1.04 it 0.16 [g]. 

In the analysis presented here, a value of the partial decay width R^ = T (W -^ cX)/r (W — > 
hadrons) is extracted from the properties of final state particles in W decays. Charm hadron iden- 
tification is based upon jet properties, lifetime information, and semileptonic decay products in 
W"*"W~ — > qqqq and W"'"W~ — > qq^£ events. The measured value of R^ is then used to deter- 
mine WcJ. 



2 Data and Monte Carlo Samples 

The data used for this analysis were recorded at LEP by the OPAL detector ^ at centre-of-mass 
energies around 183 GeV in 1997 and 189 GeV in 1998. OPALQ is a multipurpose high energy physics 
detector incorporating excellent charged and neutral particle detection and measurement capabilities. 
The integrated luminosities used for the analysis amount to 56.5 pb~ (180.2 pb~ ) at luminosity- 
weighted mean centre-of-mass energies of 182.7 GeV (188.6 GeV). In addition, 2.1 pb^^ (3.1 pb^^) of 
calibration data were collected at y/s ~ Mz in 1997 (1998) and have been used for fine tuning of the 
Monte Carlo simulation. 

To determine the selection criteria and tagging efficiencies, Monte Carlo samples were generated 
using the KoRALWl.42 Q program, which utilizes Jetset7.4 ||5[ for fragmentation, followed by a 
full simulation of the OPAL detector ^. The centre-of-mass energies were set to \/i = 183 GeV 
and 189 GeV with a W mass of Mw = 80.33 GeV. Additional samples with different centre-of-mass 
energies and W masses were used to estimate the sensitivity to these parameters. Further samples gen- 
erated with Pythia 0, Herwig [§, and EXCALIBUR [^ were used to determine the model dependence 
of the measurement. 



^The OPAL coordinate system is defined such that the 2;-axis is parallel to and in the direction of the e~ beam, the 
X-axis lies in the plane of the accelerator pointing towards the centre of the LEP ring, and the y-axis is normal to the 
plane of the accelerator and has its positive direction defined to yield a right-handed coordinate system. The azimuthal 
angle, </>, and the polar angle, 6, are the conventional spherical coordinates. 





183 GeV 


189 GeV 


Type 


qq^ 


qqqq 


qq^ 


qqqq 


Selected Events 


350 


432 


1235 


1532 


Z77 ^ qq 

qqqq 

qq^+^- 
Wez7e 


15.4 ±1.6 
0.5 ±0.6 

10.2 ±1.2 
6.2 ±1.7 


77.4 ±8.5 
12.4 ±2.8 
4.1 ±0.3 
0±1 


45.5 ±5.6 
2.5 ±0.6 

27.6 ±2.7 
23.3 ±4.1 


240.0 ± 15.9 
70.3 ±12.3 

8.6 ± 1.4 
0±1 


Total Background 


32.4 ±2.7 


93.9 ±9.0 


98.9 ±7.5 


318.9 ±20.2 



Table 1: Number of selected WW candidate events and those expected from the main background 
sources for the selection at 183 GeV and 189 GeV. The "qqqq" background source refers to processes ex- 
cluding doubly resonant W-pair production at tree level (CC03 diagrams). The entry labeled "WeI7e" 
refers to events where a singly produced W decays hadronically. The errors are the systematic errors 
on the background contributions, as evaluated in References |^l||l2[. 



The dominant background sources for this analysis are Z°/7 —>■ qq events and four-fermion final 
states which are not from W"'"W~ decays. Background samples were generated using the Pythia, 
Herwig, Excalibur, and grc4f U^ generators. More details on the treatment of the backgrounds 



can be found in References [11, 12|. 



3 Event Selection 



The WW event selections are the same as those used for the OPAL WW production cross-section 
measurements at ^/s = 183 GeV |ll| and 189 GeV |l|]. Events not selected as W+W" -^ i+uii'~Vi> 
candidates are passed to the W'^W^ -^ <^^t selection procedure and may be classified as qqel^e, 
qq/iI7^, or c^tV-j- candidates. The W"^W~ -^ qqqq selection is applied to events that fail the preceding 
selections. It has been checked that the event selection introduces a negligible bias in the flavour 
composition. In this analysis, the small fraction of selected singly produced W bosons (Wel^e) which 
decay hadronically to a charm quark is considered as signal. 

Since the identification of charm mesons relies heavily on lifetime and lepton information, the 
relevant parts of the detector were required to be operational while the data were recorded. Thus, in 



addition to the requirements which were already applied in References ||ll, 12|, it was required that 
the silicon microvertex detector and the muon chambers were fully operational. After all cuts, a total 
of 350 (1235) W"'"W~ — > qq^^ candidates and 432 (1532) W+W" -^ qqqq candidates were selected 
at 183 GeV (189 GeV). A breakdown of the different background contributions is given in Table |^. 
The expected number of background events is calculated from the background cross-sections given in 
References [^,^|; the errors correspond to the systematic errors quoted in these references. 



4 Measurement of R 



w 



After the event selection, the hadronic part (which excludes the identified lepton) of a qq^ event 
is forced into two jets using the Durham algorithm ||l^, while qqqq events are forced into four jets. 
Subsequently, a relative likelihood discriminant is calculated for each jet to separate jets originating 
from charm quarks and jets from u, d, and s quarks. This discriminant relies on jet and event shape 
properties, lifetime information, and lepton identification as described in the remaining parts of this 



section. The resulting likelihood distribution is fitted to obtain the relative fractions of uds and charm 
jets. 

In qqiPi events the two jets can be uniquely assigned to the decay of one W boson, while for qqqq 
events the four jets have to be grouped into pairs that are assumed to originate from the decay of one 
W boson. The jet pairing is chosen using a kinematic fit to find the jet combination which is most 
compatible with the production of two on-shell W bosons with a mass of 80.33 GeV. After the jet 
pairing has been performed in qq^^ and qqqq events, the jet energies are calculated in a kinematic 
fit, imposing energy-momentum conservation and equality of the masses of both W candidates. If the 
kinematic fit fails, the jet energies of each jet pair are scaled so that their sum equals the beam energy. 
The choice of the jet pairing influences the result of the present analysis only through the values of the 
jet angle 0jet in the di-jet rest frame and the di-jet angle ^w in the laboratory frame, both of which 
are calculated from the jet momenta determined by the kinematic fit. These two angles are employed 



in the likelihood functions as described in Sections 4.2 and 4.3. It has been checked in the W mass 



measurement analysis [14,15] and in an analysis concerned with properties of hadronically decaying 
W bosons ||l^ that jet properties relevant for the jet pairing algorithm and the kinematic fit are well 
modelled in the Monte Carlo simulation. 



4.1 Secondary Vertex Tag 

Weakly decaying charm hadrons have lifetimes between 0.2 ps and Ips Q, leading to typical decay 
lengths of a few hundred microns to a few millimetres at LEP2 energies. These relatively long-lived 
particles produce secondary decay vertices which are significantly displaced from the primary event 
vertex. The primary vertex in an event is found by fitting all charged tracks to a common point in 
space [|l7|. Its determination is improved by using the average beam position as a constraint with 
RMS values of 150 /zm in x, dominated by the width of the beam, about 20 /im in y, and less than 
1 cm in z. 

In the search for charm decay products, tracks in the central detector and clusters in the electro- 
magnetic calorimeter within a given jet are considered. Tracks used to find secondary vertices must 
fulfill the quality criteria used in the WW cross-section analyses [11, 12 1. Additionally, the tracks 



are required to have a momentum larger than 0.75 GeV, a signed distance of closest approach to the 
primary vertex in the rep plane (2D-impact parameter) cm < do < 0.3 cm, and an error on the 
impact parameter, ad^, smaller than 0.1cm. From the tracks which fulfill these criteria, the tracks 
which are most likely to originate from the decay products of a charm hadron are identified via an 
iterative procedure. Particles from charm hadron decays tend to have a large momentum and a large 
rapidity y = i In ( rj_ J ) (where p\\ is the particle's momentum component parallel to the jet direction 



2 ^^^ V-^~P|I 

and E is its energy) compared to fragmentation tracks, and their combined invariant mass should be 
compatible with a charm hadron mass. Therefore for each track and electromagnetic cluster which is 
not associated to any track, the rapidity y is calculated assuming the charged pion mass for tracks 
and zero mass for unassociated clusters. The track or cluster with the lowest rapidity is removed from 
the sample and the invariant mass of the remaining tracks and clusters is calculated. This procedure 
is repeated until the mass of the remaining group of tracks and clusters falls below 2.5 GeV. This 
procedure selects preferentially tracks from the decay of charm hadrons [18|. 



All selected tracks in the jet are fitted to a common vertex in three dimensions. Tracks which 
contribute more than 4 to the x^ of this fit are removed and the fit is repeated. The procedure is 
stopped when all tracks pass the x^ requirement. At least two tracks are required to remain in the fit 
for the secondary vertex finding to be successful. Once tracks are associated to a secondary vertex, the 



distance L between the primary and the secondary vertex and its error gl are computed. The decay 
length significance of a secondary vertex, L/a^^ is required to be positive (i.e. L/gl > 0), where L is 
given a positive sign if the secondary vertex is displaced from the primary vertex along the direction 
of the jet momentum. 

Jets containing a good secondary vertex are then used to tag charm decay products by means of an 
artificial neural network (ANN) |jl^ to enhance the charm purity at a given efficiency. The nine ANN 
input variables are: the decay length significance L/ai of the secondary vertex; the primary vertex 
joint probability (for a definition, see Reference [^); -E'vertex/-E'beam> where -Evertex is the sum of the 
energy of all the tracks assigned to the secondary vertex and -Ebeam is the beam energy; the distance 
of closest approach to the primary event vertex in r(j)z (3D) of the pseudo-particle formed from all 
the tracks in the secondary vertex; the normalised weighted average (weighted by the measured error) 
of all the track 2D-impact parameters within the secondary vertex; the highest momentum of the 
charged particles in the secondary vertex; the invariant mass of the vertex; the largest track 3D- 
impact parameter in the secondary vertex; and the second largest track 3D-impact parameter in the 
secondary vertex. 

The same vertex ANN is used for W^W~ -^ qqqq and W^'^W" -^ qq^£ events. The distributions 
of the two most discriminating input variables for the vertex ANN are shown in Figure |l[ Typically, 
charm hadrons have secondary vertices which are displaced from the primary vertex and charged tracks 
associated with charm particles have a lower probability to come from the primary vertex. Figure ya 
shows a comparison of the ANN output between data and Monte Carlo for ^/s = 183 — 189 GeV. The 
simulation of the OPAL tracking detectors p| requires a detailed understanding of the track-finding 
efficiency and resolution. This is achieved by tuning the Monte Carlo simulation with the calibration 
data taken at the TP resonance. Figure |3|a shows the distribution of L/a^ for secondary vertices. The 
effect of a 5 % variation of the tracking resolution used to assess the systematic uncertainty is also 
depicted. 



4.2 Lepton Tag 

About 20 % of all charm hadrons decay semileptonically and produce an electron or a muon in the 
final state. Because of the relatively large mass and hard fragmentation of the charm quark compared 
to light quarks, this lepton is expected to have a larger momentum than leptons from other sources 
(except the small contribution from semileptonic bottom decays in background events). Therefore, 
identified electrons and muons can be used as tags for charm hadrons in W boson decays. 

In a first step, the leptons are identified using ANN algorithms. The electron ANN described in 
References ||2l| , |22|| is used. Its most important input variables are the specific ionization in the central 
jet chamber (d-E/dx) of the electron candidate track and the ratio of the associated energy deposited 
in the electromagnetic calorimeter to the momentum of this track (E/p). Photon conversions are 
rejected from the sample using another dedicated ANN |2^, ^. The muon ANN |^3| is based on 
information from the tracking detectors (including dE/dx), the hadronic calorimeter, and the muon 
chambers. The lepton candidates have to be well contained in the detector (| cos 0^1 < 0.9, where 
9i is the polar angle with respect to the e~ beam direction) and are required to have momentum 
pe>2 GeV. 

The most energetic lepton candidate in each jet is used to calculate a relative likelihood discriminant 
with the goal of separating leptons in charm decays from leptons from other sources. Dedicated 
likelihood functions are set up for electrons and muons in W'''W~ — > qqqq and W'''W~ — > qqiUi 
events, using additional variables that describe properties of the lepton within a jet and enhance the 
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Tagging 
Efficiency (%) 


183 GeV 


189 GeV 


qq£l7 


qqqq 


qq^ 


qqqq 


ew^d 
ew^c 


4.3 ±0.5 

3.1 ±0.4 

3.2 ±0.4 
18.1 ±0.6 


3.8 ±0.3 

3.0 ±0.2 

3.0 ±0.2 

15.8 ±0.4 


4.9 ±0.5 

3.2 ±0.4 

3.3 ±0.4 
18.9 ±0.6 


4.4 ±0.3 

3.1 ±0.2 

3.1 ±0.2 

16.7 ±0.4 


ebgd 


6.0 ±2.7 


4.7 ±1.0 


8.0 ±2.1 


5.4 ±0.8 



Table 2: Tagging efficiencies (in %) after the WW selection from Monte Carlo simulation for u, d, s, 
and c jets in WW events and for jets from background events, for a cut on the combined likelihood 
of 0.7. The values are shown separately for events at 183 GeV and 189 GeV and for qq^I7 and qqqq 
events. The errors are statistical only. 

separation of charm jets from uds jets. The likelihood calculation is done using the method described 
in Reference |24|. The following input variables are used both for electrons and for muons: the output 
of the electron or muon ANN; the magnitude of the cosine of the polar angle of the lepton candidate 
I cos 9(\] and the magnitude of the momentum of the lepton candidate pi. 

In qqqq events, where the charge of the decaying W boson is not known, the likelihood function 
also uses the product of the lepton charge and the cosine of the W production angle, Qi cos 6'w) where 
^w is defined as the polar angle of the sum of momenta of both jets of the pair that contains the lepton. 
This quantity makes use of the strong forward-backward asymmetry of the produced W bosons. 

In qq£l7 events, where the charge of the W boson decaying leptonically is known, the likelihood 
function uses the product of the charge of the lepton tagged in the charm decay and the charge of the 
lepton from the leptonically decaying W. Finally, for electron candidates, the likelihood discriminant 
employs the output of the conversion finder ANN to suppress electrons from photon conversions. 

In Figure |2|b the likelihood output of the lepton tag is shown for data and Monte Carlo at ^/s = 
183 — 189 GeV. In this analysis precise modelling of the lepton identification is important. The lepton 
reconstruction efficiencies in the Monte Carlo simulation were tuned to match the values obtained from 
data collected at the TP resonance. Figures ^-c show the E/p and normalised dE/dx distributions 
of electron candidates. Reasonable agreement between data and Monte Carlo simulation is observed. 
The variables used for muon identification are also well described in the simulation. 



4.3 Combined Tag 



For each jet, the vertex ANN and the lepton likelihood discriminant are merged into a single quantity by 
an additional likelihood function. This combined likelihood discriminant uses two additional variables: 
the cosine of the jet angle cos 0jet in the rest frame of the W boson with respect to the W momentum 
direction and the likelihood output of the WW event selection [11, 12|. The jet angle in the di- 
jet rest frame helps to separate jets originating from up-type and down-type quarks. Due to the 
polarisation of the W bosons and the V — A nature of the weak interaction, up-type quarks tend 
to be produced preferentially in the backward direction with respect to the W flight direction, while 
down-type quarks are found more often in the forward direction. The likelihood output of the WW 
event selection suppresses false tags from background events, especially from Z'^/7 -^ bb events. 

In Figure]^ the combined likelihood output is shown, illustrating the discriminating power between 
charm and light flavoured jets. The efficiency and purity for a particular cut on the combined likelihood 



output are also depicted in Figure ^. Table |2| lists the tagging efficiencies for jets in WW and 
background events that result from a cut on the likelihood output at the optimal value of 0.7. The 
tagging efficiencies at 189 GeV are slightly higher due to the extra boost of the W bosons. These 



values are not directly used to derive R^ , but are quoted to give an impression of the performance of 
the charm tag. 



4.4 Determination of R^ 

To extract R^ from the data, a binned maximum likelihood fit is performed to the combined like- 
lihood distributions. Template histograms for the charm, light -flavoured, and non~WW background 
components are produced from the main Monte Carlo samples. They are referred to as the probability 
density functions (pdf) for the various components. Using this fitting method improves the statistical 
sensitivity of the measurement compared to simply counting the number of jets that pass a cut on the 
likelihood output. The function used to fit the data (see Figure Ha) is: 



N 



In/: = -Y^n,ln{{l-F^^^){lR^ri + [l-^,R^]Kj + F^^^n^^}, 



(1) 



1=1 



with 



i: Bin index in the combined likelihood distribution; 

Tij: Number of jets in the i bin; 

N: Number of bins in the combined likelihood distribution; 

VI: Probability for a charm jet to appear in the i bin (W ^ cX pdf); 

^uds" Probability for a uds jet to appear in the i^^ bin (W -^ light flavour pdf); 

Fbkg: Fraction of background jets; 

^bkg* Pi'obability for a non-WW jet to appear in the i^^ bin (non-WW background pdf). 



The result of the fit for i?^ is: 



and 



R. 



w 



R. 



w 



0.493 ± 0.090 at 183 GeV 



0.478 ± 0.047 at 189 GeV. 



where the errors are statistical only. The x^/(degree of freedom) of the fits are 17.2/24 at 183 GeV 
and 35.4/24 at 189 GeV. 



4.5 Systematic uncertainties 



Since the reference histograms, the efficiencies, and the background for the calculation of R^ were 
taken from a Monte Carlo simulation, the analysis depends on the proper modelling of the data 
distributions in the simulation. The sources of systematic error are listed in Table |^ and are discussed 
below. Unless otherwise noted, all changes were applied to signal and background Monte Carlo 
samples. 



Hadronisation model: The Pythia Monte Carlo sample, which uses the Lund string fragmenta- 
tion as implemented in Jetset ]5[ , was compared to a sample generated with Herwig ||8[ , which 
uses cluster fragmentation. The samples differ in their parton showering and hadronisation. 

W mass and centre-of-mass energy: As the jet pairing depends on the W mass, Mw, and the 
centre-of-mass energy, ^/s, the analysis was repeated with different W masses and at different 
centre-of-mass energies. The W mass was varied by 100 MeV, which covers both the experi- 
mental uncertainty and the differences between the current world average mass Q and the value 
used in the Monte Carlo. The centre-of-mass energy was varied by the difference between the 
luminosity-weighted centre-of-mass energies in data and the value of ^/s used in the main Monte 
Carlo samples. 

Charm fragmentation: At ^/s = Mz , the mean scaled energy of weakly decaying charm hadrons 



was found to be (xd(Mz)) = 0.484 it 0.008 [25|. No measurements of this quantity exist for W 



decays, but it is expected to be similar to that in Z decays. The changes in the charm fragmen- 
tation due to QCD effects were taken from Jetset j^. To assess the systematic uncertainty, 
(xd(Mw)) was varied by the experimental uncertainty obtained on (xniMz)), i.e. 0.008. The 



fragmentation schemes of Peterson ||22l, Collins and Spiller |27], Kartvelishvili [^ and the Lund 
group p^] were used to assess a possible discrepancy on the shape of the fragmentation function. 
The largest variation was assigned as the systematic error. 



Background cross-section: In the fit for extracting R^ , the fraction of background was varied 



according to the error on the expected fraction of background events from References [11, 12| 



Background composition: The various background sources have different probabilities to fake a 
charm jet. The background composition was varied within the uncertainties given in Table ||. 

Charm hadron fractions: The fractions of the weakly decaying charm hadrons were varied within 
their experimental errors according to Reference |25||. 



Light quark composition: The shape of the reference histogram for light-flavoured quarks de- 
pends on the relative fractions of the u, d, and s jets. The quark flavour fractions were varied 
according to the actual errors on the CKM matrix |^ so that the error covers the effect of 
assumptions concerning the CKM matrix elements on the uds template histogram. The num- 
ber of bottom hadrons from W bosons was also doubled to ensure that the measured R^ was 
insensitive to the bottom fraction in W decays. 

Vertex reconstruction: The correct modelling of the detector resolution in the Monte Carlo 
simulation is of paramount importance for the description of the probability to find a secondary 
vertex. This was checked by comparing the track parameters between calibration data collected 
at ^/s ~ Mz and Monte Carlo simulation. The sensitivity to the vertex reconstruction was 
assessed by degrading or improving the tracking resolution in the Monte Carlo. It was found 
that changing the track parameters resolution by ±5% covers the range of the observed differences 
between data and Monte Carlo (see Figure ya) in the vertex ANN input and output distributions. 

Charm hadron lifetimes: The lifetimes of the weakly decaying charm hadrons were varied within 
their experimental errors according to Reference Q. 

Charm decay multiplicity: The performance of the secondary vertex finder depends on the 
multiplicity of the final state of the charm hadrons. The relative abundance of 1, 2, 3, 4, 5, and 
6 prong charm decays was changed in the Monte Carlo within their experimental errors [P0|] . 
Similarly, the relative abundance of neutral pions was varied in the Monte Carlo within the 



experimental errors |3C]. 



Source of 


Ai?W 


Systematic Error 


183 GeV 


189 GeV 


Hadronisation Model 


0.011 


0.012 


Centre-of-mass Energy 


0.007 


0.005 


Mass of the W Boson 


0.004 


0.003 


Charm Fragmentation Function 


0.007 


0.007 


Background Cross-Section 


0.006 


0.005 


Background Composition 


0.010 


0.009 


Charm Hadron Fractions 


0.006 


0.007 


Light Quark Composition 


0.007 


0.005 


Vertex Reconstruction 


0.016 


0.017 


Charm Hadron Lifetimes 


0.003 


0.002 


Charm Decay Multiphcity 


0.010 


0.010 


Lepton Identification 


0.012 


0.014 


Lepton Energy Spectrum 


0.003 


0.003 


Branching ratio Br(c — > i) 


0.006 


0.006 


Total systematic error 


0.032 


0.032 


Statistical error 


0.090 


0.047 


Value of R^ 


0.493 


0.478 



Table 3: Summary of the experimental systematic errors, statistical errors, and central values from 
the binned maximum likelihood fit on RY^ at 183 GeV and 189 GeV. 



Lepton identification: The modelling of the input variables of the ANN for electron identification 
has been intensively studied at LEPl [^2|, while the performance of the muon tagging ANN has 
been investigated using various control samples of calibration data collected at the Z" resonance 
during each year of LEP2 |23|. The relative error of the efficiency to identify a genuine elec- 
tron (muon) has been estimated to be 4% (3%), and the relative error on the probability to 
wrongly identify a hadron as a lepton has been evaluated to be around 16 % (10 %) for electrons 
(muons) [22,^|. By studying the dE/dx distributions for electron candidates in calibration 
data taken at the Z*^ mass in 1997 and 1998, it was found that a shift of the mean dE/dx of 
0.04 keV/cm (0.3%) and a scaling by 1% of the dE/dx resolution in data was needed to correct 
a slight discrepancy between Monte Carlo and data. This effect was also taken into account 
in the assignment of the uncertainty on the electron tagging efficiency. The (mis) identification 
probabilities for electrons and muons were varied by reweighting the Monte Carlo events and 
the largest resulting change in R^ was taken as the systematic error. It was verified that the 
variation of efficiencies and dE/dx resolution was sufficient to describe the possible differences 
between data and Monte Carlo (see Figures |^-c) . 

Lepton energy spectrum from semileptonic charm decays: The energy spectrum of the lepton 
in the rest frame of the weakly decaying charm hadron was reweighted from the spectrum 
implemented in Jetset to the ISGW Q and ACCMM |3^ models; the parameters pf and nis 
of the ACCMM model were varied in the range recommended in Reference 12^ 



Branching ratio Br(c -^ i): The branching fraction for inclusive charm decays to leptons was 
varied within the range Br(c -^ i) = 0.098 it 0.005 [25| by reweighting the Monte Carlo events. 



The size of the systematic errors is very similar for the samples at 183 GeV and 189 GeV. The dominant 
systematic errors are those associated with the vertex reconstruction and lepton identification. 
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CKM matrix element 


Value 


Kd 

\Vus\ 

\Vnk\/\Vcb\ 

l^cdl 

iKbl 


0.9735 ± 0.0008 
0.2196 ±0.0023 
0.090 ±0.025 
0.224 ±0.016 
0.0402 ±0.0019 



Table 4: Values of the CKM matrix elements used in the calculation of |V^s| from R^ . The values 
are taken from Reference |^]. 

5 Results 



The result of the charm tag analysis is 



R^ = 0.493 ± 0.090 (stat.) ± 0.032 (syst.) at 183 GeV 



and 



R^ = 0.478 ± 0.047 (stat.) ± 0.032 (syst.) at 189 GeV. 



The combined value of i?^ obtained at ^/s = 183 GeV and 189 GeV is 



R^ = 0.481 ± 0.042 (stat.) ± 0.032 (syst.), 
where the systematic uncertainties are treated as being fully correlated. 



6 Extraction of IV^ 



cs 



It is possible to extract |T4s| from R^ using the relation 



R. 



\Vnd\^ + iKsP + iKbP ± ll^cdP ± iKsP ± l^cbP' 



which leads to 



Wr, 



i?w 



lKdP±|KsP±|KbP) 



iKdP - |l/cbP. 



(2) 



(3) 



Using the values of the other CKM matrix elements listed in Table |^, one gets 

iKsl = 0.93 ± 0.08 (stat.) ± 0.06 (syst.) ± 0.004 (CKM), 
where the last error is due to the uncertainty on the CKM matrix elements used in the calculation. 

A more precise value for |I4s| can be obtained based on the Standard Model prediction of the 
decay width of the W boson to final states containing charm quarks: 



r(w^cX) 



^0|^(|KdP±|Ks|^±|KbP), 



(4) 



with 



6V27r 



(705 ± 4) MeV @ and the color factor C given by @ 



^ ^ 3 , ^ ^ MMw) + i.409^ii^ 



12.77 



"s (^w: 



IT 



■K^ 



^o 



(5) 
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r (W — > c X) can be evaluated from R^ using the relation 

r (W ^ c X) = i?^ Br ( W ^ hadrons) T^^ . (6) 

Using the PDG values Br (W -^ hadrons) = 0.6848 ± 0.0059 and0 T^t = (2-12 ± 0.05) GeV |] 
together with the measurement of R^ presented in this letter, one gets 

r (W ^ cX) = [698 ± 61 (stat.) ± 46 (syst.) ± 18 (ext.)] MeV, 

where the last error is due to the uncertainties on the external parameters F^^ and Br ( W -^ hadrons) . 
This corresponds to 

iV^dl^ + iV^csl^ + iV^cbl^ = 0.990 ± 0.087 (stat.) ± 0.065 (syst.) ± 0.026 (ext.). 

Subtracting the off-diagonal CKM matrix elements, using the values listed in Table ^, results in 

l^csl = 0.969 ± 0.045 (stat.) ± 0.034 (syst.) ± 0.013 (ext.) ± 0.004 (CKM), 

where the last error is due to the uncertainties on the CKM matrix elements used in the calculation. 
The determination of |V^s| from Equations (0) and IM) only involves the CKM matrix elements for 
the transition W — > c X and well measured parameters such as Pj^^ and Br ( W -^ hadrons) . It is 
therefore more precise than the method which relies on Equation (^). Both results of |V^s| presented 
here are compatible with the direct determination from D semielectronic decays which yields |V^s| = 

1.04 ±0.16 a. 



7 Summary 

Using data taken with the OPAL detector at LEP at centre-of-mass energies of 183 GeV and 189 GeV, 
the ratio R^ = P (W —>■ cX)/P (W — > hadrons) has been measured using a technique based on jet 
properties, lifetime information, and leptons from charm decays. The combined result is 

R^ = 0.481 ± 0.042 (stat.) ± 0.032 (syst.), 

in agreement with the Standard Model expectation of 0.5. This result is consistent with the Standard 
Model prediction that the W boson couples to up and charm quarks with equal strength and it is in 



agreement with the recent measurements of DELPHI |34| and ALEPH 



Using the results from direct measurements of the other CKM matrix elements, |V^s| can be 
calculated from R^^ , yielding |V^s| = 0.93 ± 0.10. A more precise value for iVcsl can be derived from 
the measurements of R^ , Br (W -^ hadrons) , and P^^^ : 

\Vcs\ = 0.969 ± 0.058. 

This can be compared to a more indirect determination of |V^s| by OPAL from the ratio of the W 
decay widths to qq and iPe, which yields jFcsl = 1-014 ± 0.029 (stat.) ± 0.014 (syst.) H under the 
assumption that the W boson couples with equal strength to quarks and leptons. It should be noted 
that the ratio Br (W — > qq)/Br (W — > iVi) is not used in the analysis presented here, so that these are 
independent determinations of |V^s|- 

^ The most precise measurements of Fj^t are determined from the quantity 7?. — ""(pp-g^ + y ^JoJ^T-) ■ In |33| , 
the coupling of the W boson to Ught quarks is needed as input; however, direct measurements of the corresponding CKM 
matrix elements are so precise that it introduces a negligible uncertainty on Fj^t- 
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(b) 183-189 GeV 
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Figure 1: Distributions of the most sensitive input variables of the vertex ANN at ^/s = 183—189 GeV: 
(a) decay length significance L/a^ and (b) primary vertex joint probability. Points with error bars 
represent the data (183 GeV and 189 GeV samples combined). Histograms represent the Monte Carlo 
expectation. 
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Figure 2: Distributions of (a) the vertex ANN output and (b) the hkehhood output of the lepton 
tag at -y/i = 183 — 189 GeV. Points with error bars represent the data (183 GeV and 189 GeV samples 



combined). Histograms show the results of the fits to the combined likelihood described in Section i.4. 
If no vertex or lepton tag information is available, the corresponding variable is assigned the value of 
zero. 
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Figure 3: Comparison of discriminating variables between calibration data taken at the Z^ resonance 
and Monte Carlo prediction. Points with error bars represent data from the Z^ calibration runs. 
Histograms show the expectation from the Monte Carole simulation, (a) Decay length significance 
L/aL for data and Monte Carlo. The dotted line shows the effect of the ±5 % variation of the tracking 
resolution used to assess the systematic uncertainty on the vertex tag. (b) E/p and (c) normalised 
dE/dx for electron candidates from data and Monte Carlo. In these plots, electrons must satisfy the 
pre-selection described in the text and have an electron ANN output greater than 0.5. The dotted 
lines show how the variation of the detection (mis)identification and resolution for electron candidates 
affects the E/p and dE/dx distributions. 
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Figure 4: (a) Output of the combined likelihood used to tag charm hadrons at 183 — 189 GeV. Points 
with error bars represent the data (183 GeV and 189 GeV samples combined). Histograms represent 
the results of the fits described in Section [4.4| . (b) Efficiency and purity computed after the WW 
selection as a function of a cut on the combined likelihood output. 
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